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Abstract 

A measurement of the A^ lifetime using the decay A^ — > ]/tpA in proton-proton col- 
lisions at v^ = 7 TeV is presented. The data set, corresponding to an integrated lumi- 
nosity of about 5 fb^ , was recorded with the CMS experiment at the Large Hadron 
Collider using triggers that selected dimuon events in the ]/xp mass region. The A^ 
lifetime is measured to be 1.503 ± 0.052 (stat.) ± 0.031 (syst.) ps. 
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1 Introduction 

The study of b baryons is a necessary ingredient to understand b-hadron phenomenology. 
The heavy-quark expansion model of nonperturbative quantum chromodynamics provides a 
framework for predicting properties of heavy-flavour hadrons, including their lifetimes. Here, 
a simple description is used where the heavy b quark is surrounded by a light quark or di- 
quark system. Estimates can be made of the lifetime and of the ratio of lifetimes between par- 
ticles sharing the same heavy-quark flavour JlHSl. The early calculation predicted a spread 
of the lifetimes of order 5% among all b hadrons [6J and the ratio of the A^ and B*^ life- 
times, t^o/tbo, to be greater than 0.90 tZJ- The initial measurements of the b-baryon lifetime 
were generally lower than predicted [Zl, w^hich motivated more refined calculations. This re- 
sulted in predictions as low as Tao/tbo = 0.86 ± 0.05 ||8|. An overview of the current state 

b 

of the predictions and measurements can be found in Ref. IZl. Measurements of the A^ life- 
time prior to 2011 can be found in Refs. Il9]- t25l . More recent measurements of the A^ lifetime 
include: 1.537 ± 0.045 ± 0.014 ps from CDF 1261, 1.303 ± 0.075 ± 0.035 ps from DO [27J, and 
1.449 ± 0.036 ± 0.017 ps from ATLAS ||28|, where the first uncertainties are statistical and the 
second are systematic. 

In this paper, a measurement of t^o is presented, using the decay A^ -^ ]/^A, with A — > p/r^ 

etnd]/xp —7- }i^H^- The kinematically analogous channel B*^ — )■ J/i/^Kg, with Kg — ?► tt+tt^, is used 
as a cross-check, with selection criteria similar to those for the A^ analysis. Charge-conjugate 
states are assumed throughout this paper. The measurement is made using proton-proton col- 
lision data at -^s = 7TeV recorded by the Compact Muon Solenoid (CMS) experiment oper- 
ating at the Large Hadron Collider (LHC). The data set for this measurement was collected in 
2011 using J/t/7-enriched dimuon triggers, and corresponds to an integrated luminosity of about 

5fb"^ mi. 



2 CMS detector 

The central feature of the CMS apparatus is a superconducting solenoid of 6 m internal diam- 
eter. The main subdetectors used for the analysis are the silicon tracker, consisting of silicon 
pixel and silicon strip layers, and the muon system. The tracker is immersed in a 3.8 T axial 
magnetic field of the superconducting solenoid. The pixel tracker consists of three barrel layers 
and two endcap disks at each barrel end. The strip tracker has 10 barrel layers and 12 endcap 
disks at each barrel end. The tracker provides an impact parameter resolution of ~15 jim and 
a transverse momentum (px) resolution of about 1.5% for 100 GeV particles. Charged hadrons, 
including pions and protons, are not explicitly identified by their type. Muons are measured 
in gas-ionisation detectors that are embedded in the steel return yoke outside the solenoid. In 
the barrel, there is a drift tube system interspersed with resistive plate chambers, and in the 
endcaps there is a cathode strip chamber system, also interspersed with resistive plate cham- 
bers. The first-level trigger used in this analysis is based on the muon system alone, w^hile the 
high-level trigger uses additional information from the tracker. A detailed description of the 
CMS detector can be found in Ref. [,30 J . 

The CMS experiment uses a right-handed coordinate system, with the origin at the nominal 
interaction point, the x axis pointing towards the centre of the LHC ring, the y axis pointing 
up (perpendicular to the plane of the LHC ring), and the z axis along the anticlockw^ise-beam 
direction. The polar angle 9 is measured from the positive z axis and the pseudorapidity is 
defined hy rj = — ln[tan(0/2)]. The azimuthal angle is measured from the positive x axis in the 
plane perpendicular to the beam. 



2 3 Event selection and efficiency modelling 

3 Event selection and efficiency modelling 

Dimuon triggers optimised for selecting events with ]/xp candidates are used. The trigger re- 
quires two oppositely charged muons with an invariant mass compatible with the ]/xp-meson 
mass. The dimuon candidate must also be found in the central region (dimuon rapidity ly^+^f- 1 < 
1.25), which is the region with the best impact parameter and dimuon invariant-mass resolu- 
tion. With increasing instantaneous luminosity, the trigger requirements were adjusted several 
times during the data-taking period. In the course of data taking, the dimuon mass w^indow^ 
was changed from 2.5^.0 GeV to 2.30-3.35 GeV, while the minimum transverse momentum of 
the dimuon candidate was increased from 6.5 GeV to 13 GeV. Additional requirements were 
also added during this period: the distance of closest approach between the two muons was 
required to be less than 0.5 cm, the vertex-fit x^ probability of the two muons was required to 
be greater than 0.5%, and the two muons were required to bend away from each other in the 
tracker. 

The four charged particles (j-i^ii^prc^) in the decay channel A^ — > ]/ipA allow for a full re- 
construction of the A^ baryon. The selection requirements are chosen to maximise the ratio of 
signal yield to the square root of the signal-plus-background yield. Events w^ith two oppositely 
charged muons are selected. The muons are reconstructed using information from the tracker 
and the muon system. The muon candidates must be within the kinematic acceptance of the 
detector by demanding the muon transverse momentum pj satisfies pj > 3.5 GeV for muon 
pseudorapidity \t]^\ < 1.2, and pj > 3.0-3.5 GeV for 1.2 < \f]^\ < 1.6, where the pj threshold 
decreases linearly as a function of |f/'' |. The muon candidates are also required to have a track 
X^ per degree of freedom less than 1.8, at least 11 tracker hits, at least two hits in the pixel sys- 
tem, a match to at least one track segment in the muon system, and a transverse (longitudinal) 
impact parameter less than 3 cm (30 cm) with respect to the primary event vertex. The prelim- 
inary choice for the primary vertex is the one with the highest sum of the squares of pj of the 
tracks associated with it. The two muons are then used to form a ]/xp candidate. 

A A candidate is formed using oppositely charged tracks with the proton candidate required to 
have pT>1.8GeV and the pion candidate px>0.46GeV. The higher-momentum track is taken 
as the proton. 

Kinematic vertex fits are used to identify the A^ candidates ||3T1 . An initial unconstrained fit 
is used to measure three selection parameters. First, the proton-pion invariant mass (m-pjj^-) is 
found from the A candidate tracks and is required to be within 6 MeV of the world-average A 
mass [32J. The A candidates are rejected if the dipion invariant mass nij^+j^- is within 10 MeV 
of the Kg mass II32I , when the proton candidate is assigned a pion mass. Then, the dimuon 
invariant mass (m„+„-) is determined from an unconstrained fit to the J/xp candidate tracks 
and is required to be within 300 MeV of its world-average value |32|. Finally, a pointing angle 
is measured for the A candidate, w^hich is defined as the angle between the momentum of 
the reconstructed A particle and the vector between its production and decay vertices. The 
pointing angle is required to be less than 0.015 rad. 

To determine the proper decay time f of a A^ candidate, another kinematic vertex fit is per- 
formed where the two muon tracks and the A candidate are constrained to come from a com- 
mon vertex, with the A and ]/xp candidate masses constrained to their respective nominal val- 
ues p2| . The A vertex-fit x^ probability must be greater than 10%. The separation between the 
A vertex and the eventual primary vertex (as defined below) must be larger than lOc^, where cr 
is the calculated uncertainty in the relative position. In addition, the A vertex must be at least 

3 cm away from the mean pp collision position in the transverse plane. 



Multiple pp interactions per bunch crossing are present in the data. The reconstructed event 
vertex with the smallest impact parameter in the z direction to the A^ candidate trajectory is 
selected as the primary production vertex. The position of the primary vertex is recalculated 
excluding the tracks from the A^ decay if they were used for the initial primary-vertex recon- 
struction. The proper decay time t is found for each A^ candidate by calculating the ratio of the 
decay length and the momentum of the A^, divided by the world-average A^ mass II32II . 

As a cross-check, the B^ lifetime, Tgo, is determined using the B*^ — > J/i/?Kg channel. A com- 
pletely analogous procedure is followed to find the B*^ candidates, replacing the proton-track 
hypothesis with a pion-track hypothesis and fitting for a Kg candidate instead of a A candidate. 
The Kg candidates are formed assuming both tracks are pions (the higher-momentum one is re- 
quired to have pT>l-8 GeV, the lower-momentum one to have pT>0.5 GeV), with m„+jj-- within 
12 MeV of the world-average value [32] . A candidate is vetoed if replacing the pion-mass hy- 
pothesis by the proton mass yields an invariant mass ntpj^-- within 10 MeV of the A mass [|32| . 
The pointing angle and other kinematic requirements are the same as those for the A^ selection. 

Simulated event samples are used to model the signal and background distributions. The 
PYTHIA 6.422 1331 event generator is used, with the A^ lifetime fixed to 1.425 ps 1321 and the b- 
hadron decays described by the evtgen 9.1 simulation package [341. The particle propagation 
and detector simulation is performed using Geant4 9.4 [35 1, and the events are fully recon- 
structed with the same software as the data. The differences between the reconstructed and 
generated values for the proper decay time, the flight length, and the candidate A^ momentum 
are found to be compatible with zero in the simulated sample. The sideband-subtracted data 
distributions in each of the selection variables are compared to those in simulated data and 
found to be consistent within their uncertainties. The overall reconstruction and selection effi- 
ciency as a function of the proper decay time t is determined from the simulated signal samples 
by calculating the ratio of the numbers of reconstructed and generated A^ (B*^) candidates in 
bins of proper decay time. Figure [l] displays the reconstruction and selection efficiencies mea- 
sured from simulation as a function of the proper decay time for the B*^ (top) and A^ (bottom). 
The efficiencies are consistent with being independent of the proper decay time, as shown by 
the solid horizontal line in each plot. 

4 Proper decay time fit 

Unbinned extended maximum-likelihood fits are performed to determine the A^ and B° life- 
times. The input variables are the invariant mass m, proper decay time t, and its uncertainty at, 
calculated per candidate from the kinematic vertex fit using full error propagation. The likeli- 
hood fit is implemented using the ROOFiT 3.53 package II36I . The likelihood function (ignoring 
the normalisation terms for simplicity) is 

■^ = n ^sig ■ G2(m;; msig, (7,„j, (Tm^,/) ■ e"*''"'^'8 (g) G{ti) ji, S ■ at,i) 
i 

+ Nprompt ■ Pirn; a) ■ G{iu V, S ■ at^i) (1) 

+ Nnonprompt • P{mi) a) ■ e-'/"""-?"-?' ^ G(i,-; ^, S ■ atj) , 

where the index i goes over the events, N^ig is the number of signal events, Nprompt is the 
number of prompt background events not coming from b-hadron decays, and Nnonprompt is 
the corresponding number of nonprompt background events coming from b-hadron decays. 
The prompt background is dominated by ]/ip mesons directly produced in the pp collision, 
while the nonprompt background is dominated by ]/xp mesons from decays of b hadrons. In 
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Figure 1: The overall reconstruction and selection efficiency (in arbitrary units) as a function of 
proper decay time for B*^ (upper two plots) and A^ (lower two plots). The second and fourth 
plots shows the B*^ and A^ efficiencies, respectively, for smaller proper decay times f < 1 ps. 
The horizontal solid lines show the average efficiency, while the dashed lines display the result 
of the fits to the function given by Eq. (|2]|. 



both cases, the ]/xp mesons are combined with real or misidentified A candidates from the 
event. The parameters Tgig and Tnonprompt denote the lifetime of the signal and of nonprompt 
background, respectively. In Eq. ([l}, the G2 function is the sum of two Gaussians with a 
common mean mgig and widths o;„i and a,„2, and the parameter / denotes the relative fraction of 
the area of the two Gaussians. The function G refers to a single Gaussian describing the detector 
lifetime resolution, w^ith a mean fi and a width S ■ dfj, where S is a scale factor determined 
from the fit. This resolution function is common to all three likelihood components since the 
Ct distributions do not differ significantly. The effect of using this simplifying assumption on 
Ct is evaluated as a systematic uncertainty and found to be negligible. For the background 
components, the invariant-mass m distribution is parameterised by a normalised first-degree 
polynomial of slope a, P{m;a). The prompt and nonprompt backgrounds share the same slope. 
The maximum-likelihood fit to the data is performed allowing all parameters to vary. 

5 Results 

Projections of the invariant-mass and proper decay time distributions and the results of the fits 
for the B*^ and A^ are shown in the upper panels of Figs.0andpl respectively. The lower pan- 
els in each figure give the proper decay time projections and fit results for low-mass sideband 
(left), signal (center), and high-mass sideband (right) regions of the invariant-mass distribution. 
The signal region is defined to be within 2cr of the mass peak, where a is the mass resolution 
obtained by integrating the double-Gaussian signal function with its parameter values deter- 
mined from the fit. The low-mass sideband region goes from 5.15 (5.4) GeV to within 3cr below 
the peak, and the high-mass sideband region runs from 3a above the peak to 5.75 (6.0) GeV for 
the B^ (A^). The lower plot in each panel of Figs. and p^ displays the pull distribution for the 
corresponding data and fit results shown in the upper plot. 

The results from the fits are summarised in Table [ij w^here the uncertainties are statistical only. 
The measured B'' lifetime shown in the table is consistent with the w^orld-average value of 
1.519 ± 0.007 ps ||32| . The mass values for both A^ and B'^ given in the table compare w^ell with 
the world-average values [32|. 

The fitting procedure is validated by studying simulated pseudo-experiments in which the 
proper decay time distributions are generated using different lifetime values. The resulting 
lifetime measurements found from the fit are compatible with being unbiased, and the width 
of the pull distribution, (measured value — input value) /uncertainty, is consistent with 1.0. 

Table 1: Summary of the fit results for A^ and B*^ w^ith their statistical uncertainties. 



Hadron Ngig m (MeV) t (ps) 

"Ag 1013 ±40 5619.7 ±0.5 1.503 ±0.052 

B'' 6772 ±87 5278.9 ± 0.2 1.526 ±0.019 



Sources of systematic uncertainty are detector alignment, efficiency as a function of proper 
decay time, event selection, and fit model. To estimate possible effects due to uncertainties in 
the alignment, nine different simulated samples with distorted geometries are produced and 
analysed f37]. The lifetime difference between the nominal result and the sample that produces 
the largest deviation (scaled to the estimated residual misalignment present in the detector) is 
taken as the systematic uncertainty from this source. 

Since the overall efficiency, determined through simulation, is consistent w^ith being indepen- 
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Figure 2: Projections of the invariant-mass and proper decay time distributions and the results 
of the fit are shown for the B*^ decay in the upper panels. The dark solid lines give the results 
of the overall fit to the data. The lighter solid lines are the signal contributions, and the dashed 
and dotted lines show the prompt and nonprompt background contributions, respectively. The 
lower panels display the proper decay time projections for the low-mass sideband (LSB, left), 
the signal (SR), and the high-mass sideband (HSB, right) regions defined in the text. The lower 
plots in each panel give the corresponding pull distributions for the data and fit results shown. 
All plots are from the same fit. 
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Figure 3: Projections of the invariant-mass and proper decay time distributions and the results 
of the fit are shown for the A^ decay in the upper panels. The dark solid lines give the results 
of the overall fit to the data. The lighter solid lines are the signal contributions, and the dashed 
and dotted lines show the prompt and nonprompt background contributions, respectively. The 
lower panels display the proper decay time projections for the low-mass sideband (LSB, left), 
the signal (SR), and the high-mass sideband (HSB, right) regions defined in the text. The lower 
plots in each panel give the corresponding pull distributions for the data and fit results shown. 
All plots are from the same fit. 



8 6 Summary 

dent of the proper decay time, as shown in Fig. Ilj no efficiency correction is used in the life- 
time result. Nevertheless, the effect of a possible proper-decay-time-dependent efficiency is 
included as a systematic uncertainty. To this end, the efficiency is included in the likelihood 
function in Eq. ([T]) as a function of the measured proper decay time. The difference between the 
lifetime found using a constant efficiency (i.e. no efficiency in the likelihood) and that found 
using the fitted efficiency function is taken as a systematic uncertainty. The efficiency is fit to 
the function 

e{t) = po-{l + p,t+:^^!^y (2) 

where the free parameters p, are determined from the fit. The results of fitting this function to 
the B^ and A^ efficiencies are shown by the dashed lines in Fig. \\\ 

To account for possible biases from the event selection criteria, a systematic uncertainty is cal- 
culated from the difference betw^een the observed and expected values in simulated events 
using the full analysis chain. 

Simulated pseudo-experiments produced with different input parameter values and different 
modelling of the fit functions are used to estimate the corresponding systematic uncertainties 
on the lifetime measurement. Variations include: using two different lifetimes to describe the 
nonprompt background to control the possible presence of a second nonprompt background 
component, varying the prompt and nonprompt background, varying the lifetime of the non- 
prompt background to control possible correlation to the lifetime of the background, and using 
larger per-event uncertainties for the background components to control a possible mismod- 
elling of the resolution. Extending the likelihood function in Eq. ^ from a single detector life- 
time resolution at for all three components to individual resolutions per component showed a 
negligible effect. Table [2]gives the systematic uncertainties from the four sources and their sum 
in quadrature, w^hich is taken as the overall systematic uncertainty on the A^ lifetime measure- 
ment. 

Table 2: Summary of the systematic uncertainties in the A^ lifetime measurement. 



Source Systematic uncertainty (ps) 



Alignment 


0.005 


Efficiency 


0.030 


Event selection 


0.005 


Fit model 


0.004 



Total 0.031 



Checks are also performed to see if any detector effect leads to a systematic deviation not cov- 
ered by the ones previously discussed. This is done by dividing the data sample into parts, 
where each of the partitions is expected to give the same results. The checks are performed 
with azimuthal angle, pseudorapidity, transverse momentum, run era, muons bending away 
or towards each other, and number of primary interaction vertices. No statistically significant 
effects are seen. For the A^ channel, we also remove the Kg-mass-veto requirement on the A 
candidate and find a negligible effect on the lifetime result. 

6 Summary 

A measurement of the A^ lifetime has been presented using the decay A^— >J/i/?A in pp col- 
lisions at a/s = 7TeV with the CMS detector. From a data set corresponding to an inte- 



grated luminosity of about 5fb ^, the A° lifetime is found to be t^o = 1.503 ± 0.052 (stat.) ± 

0.031 (syst.)ps. The kinematically similar decay B*^— >J/t/7Kg was used as a cross-check, confirm- 
ing that no efficiency correction was needed. The A^ lifetime result is in agreement with the 
world-average value of 1.425 ± 0.032 ps |32| and has a precision comparable to that of other 
recent measurements II26H28I . As illustrated in Fig. Ill this new^ result confirms the tendency 
of the more recent measurements that give larger lifetimes, in better agreement with the early 
theoretical predictions [|6l[Z|. 
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Figure 4: Evolution over time of the A^ lifetime measurements (left scale) |9]- l28li and the ratio 
of the measurements to the 2012 world-average B*^ lifetime Tgo [32J (right scale). The values 
shown as open circles were not included in the Particle Data Group (PDG) 2012 average II32I , 
displayed as the grey band, w^hile those shown as filled circles were included. The result of this 
analysis is shown by the open square. The inner error bars represent the statistical uncertain- 
ties, and the outer error bars show the combined statistical and systematic uncertainties added 
in quadrature. Where needed, points have been shifted slightly along the time axis to enhance 
clarity. 
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